A cosmological solution of Regge calculus 



(N 



13 



Adrian P Gentle 

Department of Mathematics, University of Southern Indiana, Evansville, IN 47712 
E-mail: apgentle@usi.edu 



o 

(N 

£jq| Abstract. We revisit the Regge calculus model of the Kasner cosmology first 

^ . considered by S. Lewis. One of the most highly symmetric applications of lattice 

gravity in the literature, Lewis' discrete model closely matched the degrees of freedom 
of the Kasner cosmology. As such, it was surprising that Lewis was unable to obtain 
the full set of Kasner-Einstein equations in the continuum limit. Indeed, an averaging 
procedure was required to ensure that the lattice equations were even consistent with 
Q" 1 the exact solution in this limit. We correct Lewis' calculations and show that the 

5^ | resulting Regge model converges quickly to the full set of Kasner-Einstein equations in 

toO the limit of very fine discretization. Numerical solutions to the discrete and continuous- 

time lattice equations are also considered. 
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C\ 1. Introduction 

The discrete formulation of gravity proposed by T. Regge in 1961 pQ has been deployed 
^ ; in a wide variety of settings, from probing the foundations of gravity and the quantum 

realm [2} [3] to numerical studies of classical gravitating systems [21 II]- Regge calculus 
continues to be used in new and diverse ways; recent examples include Ricci flow [5] 
and as an explanation of dark energy [6]. 

In this paper we re-examine the Regge calculus model of the vacuum Kasner 
cosmology first considered by Lewis [TJ, with the goal of gaining insight into the 
continuum limit of this discrete approach to gravity. In a general setting the structural 
differences between a continuous manifold and a discrete simplicial lattice lead to 
difficulties in directly comparing the Regge and Einstein equations or their solutions, 
with a single Regge equation per edge in the lattice compared with ten Einstein equations 
per event in spacetime. We expect many more simplicial equations than Einstein 
equations in a general simulation, and some form of averaging must be expected before 
the equations (or their solutions) can be compared. 

Lewis [7] studied both the Kasner and spatially flat Friedmann-Lemaitre- 
Robert son- Walker (FLRW) cosmologies using a regular hypercubic lattice. We only 
consider the Kasner solution in this paper, where the high degree of symmetry, without 
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the added complication of matter, allows explicit examination of the Regge equations 
in the continuum limit. By aligning the degrees of freedom of the lattice with the 
continuum metric components, Lewis was able to avoid the issue of averaging and make 
direct comparisons between the Regge equations and the Kasner-Einstein equations in 
the continuum limit. Unfortunately, Lewis was only able to recover one of the four 
Einstein equations in this limit, and even this was only possible after the equations were 
carefully averaged [7J. Without this averaging it is not clear that the equations obtained 
by Lewis actually represent the Kasner cosmology. 

We show that Lewis neglected a vital portion of the simplicial curvature arising 
from the two-dimensional spacelike lattice faces that lie on constant time hypersurfaces. 
The Kasner cosmology has zero intrinsic curvature on constant time hypersurfaces, 
so the lattice curvature concentrated on spacelike faces - measured on a plane with 
signature "-+" orthogonal to the face - make an important contribution to the total 
lattice curvature. We show that when these curvature terms are included, the discrete 
equations exactly reproduce the Kasner-Einstein equations in the limit of very fine 
triangulations without the need to average. 

In addition to reconsidering Lewis' analytic work on the spatially flat, anisotropic 
Kasner cosmology [7J, we construct numerical solutions to the discrete and continuous- 
time lattice equations. This builds on the previous work of Collins and Williams [8] and 
Brewin [9] on highly symmetric, continuous time, closed FLRW cosmologies, and the 
(3+l)-dimensional numerical study of the Kasner cosmology by Gentle [10] with discrete 
time and coarse spatial resolution. We begin by briefly describing the continuum Kasner 
solution. 

2. The Kasner cosmology 

The Kasner solution [TT] is a vacuum, homogeneous, anisotropic cosmological solution of 
the Einstein equations with topology R x T 3 . Appropriate slicing creates flat spacelike 
hypersurfaces, while the global topology allows non-trivial vacuum solutions of the 
Einstein equations. 

The Kasner metric may be written in the form [T2] 



where the functions /, g and h are determined by the vacuum Einstein equations 



ds 2 = -dt 2 + f{t) 2 dx 2 + g{t) 2 dy 2 + h{t) 2 dz 2 
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Figure 1. The section of a world-tube joining a rectangular prism to its future 
counterpart. Homogeneity implies that an observer will fall freely along the centre 
of the worldtube, providing a convenient coordinate system from which to view the 
lattice. 



Note that the equation G u = is a first integral of the remaining equations. The Kasner 
metric components are 

f{t) = t 2 v\ g{t) = t 2 ^ 1 h(t) = t 2 ™, 

where the Kasner exponents pi are unknown constants. With this choice of metric 
functions the vacuum Kasner-Einstein equations reduce to two algebraic constraints, 

p\ +pI +pI = pi +p 2 +P3 = i, 

leaving a one parameter family of Kasner cosmologies. 

The Kasner solutions are the basis of the Mixmaster cosmologies [13], which may 
be regarded as a series of Kasner-like epochs undergoing an infinite series of "bounces" 
from one set of Kasner exponents to the next. It is conjectured that these asymptotic 
velocity term dominated models embody the generic approach to singularity in crunch 
cosmologies, and it has been shown that the bounces represent a chaotic map on the 
Kasner exponents [HI [15] . 



3. A homogeneous, anisotropic spacetime lattice 

We follow Lewis and build a discrete approximation of the Kasner spacetime using a 
highly symmetric lattice of rectangular prisms. The regularity of the lattice implements 
homogeneity, while the rectangular prisms allow a degree of anisotropy. The complete 
four-geometry is constructed by extruding the initial three-geometry forward in time 
and filling the interior with four-dimensional prisms. 

Each flat T 3 hypersurface consists of a single rectangular prism with volume Xit/iZi, 
where opposing faces are identified to give the global topology. This prism is subdivided 
into n 3 regular prisms with edge lengths u^Vi and Wi, with 

Ui Zi , . 

Ui = — , Vi = — , Wi = — , (5) 
n n n 

and where the subscript i labels the Cauchy surface at time ij. 
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The three-geometry is joined to a similar structure at time t = = tj + At i; where 
the prisms have edge lengths u i+ i, v i+ i and w i+ i. This structure is shown in figured! 
Time evolution of the initial surface maintains homogeneity, so within the worldtube of 
each prism there exists a local freely-falling inertial frame. In the coordinates of this 
frame the coordinates of vertex A in figure [T] can be written as 

ft Y± 

V*' 2 ' 2 ' 2 J 

and similarly, the coordinates for vertex A + (the counterpart of A on the next 
hypersurface) are 

2 ' 2 ' 2 )' 

The spacetime interval along the time-like edge joining A and A + is defined to be mf < 0, 
and thus 

m? = -At, 2 + \ (An} + Av} + Aw}) , (6) 

where the difference operator is defined as AU = — k. Note that the requirement 
that m 2 < implies a restriction on the choice of At, for a given value of the "resolution" 
parameter n. Homogeneity guarantees that identical expressions hold for all timelike 
cd ges joining the spacelike hypersurfaces labeled t{ and t^\. 

The discrete spacetime curvature in Regge calculus is manifest on the two- 
dimensional faces on which three-dimensional blocks hinge [1] , and is represented by the 
angle deficit (the difference from the flat space value) measured in the plane orthogonal 
to the face. There are two distinct classes of two-dimensional faces, or hinges, in 
the lattice: timelike areas formed by evolving a spacelike edge forward to the next 
hypersurface, and the rectangular faces of the prisms on t = constant slices. 

The timelike trapezoids formed when the spatial edge Ui is carried forward in time 
from the hypersurface labeled U to ti+i, face ABA + B + in figure (U has area 

Af = - (u i+1 + m) ^AvZ-Aml 

which is real since m 2 < 0. Likewise, the spacelike hinge ABCD shown in figure (TJ 
consisting of the edges U{ and Uj, has area 

Af = mvi, 

with the other spacelike and timelike hinges defined similarly. 

Turning now to the curvature about these faces, we note that there are four distinct 
three-dimensional prisms which hinge on the timelike face ABA + B + . Each of these is 
formed by dragging one of the two-dimensional faces containing the edge AB forward 
in time. This includes the prisms ABCDA + B + C + D + and ABEFA + B + E + F + . The 
remaining two prisms hinging on ABA + B + are not shown in figure [H 

The homogeneity of the lattice ensures that the four hyper-dihedral angles that 
surround the hinge ABA + B + are the same. Denoting each of these angles as 9f f , the 
angle defect (or deficit angle) about the timelike face ABA + B + is 

ef = 2tt - 49 f , 
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which measures the deviation of the total angle from the flat-space value of 2ir. The 
hyper-dihedral angle Of 1 is 



cos 9. 



xt 



v/(4Atf - Avf) (4A*f - Awf) 

with analogous definitions for Of* and 9f, the hyper-dihedral angles about the remaining 
classes of spacelike hinge. 

To measure the deficit angles about the spacelike faces, consider the curvature e® y 
about the hinge ABCD. Since this hinge is spacelike, the hyper-dihedral angles are 
boosts in the plane with signature — h orthogonal to the hinge. The deficit angle for a 
spacelike hinge is [IB] 

k 

where the summation is over all boosts 0& which surround the hinge. The boost between 
the two three-dimensional prisms which hinge on ABCD, namely ABC D A + B + C + D + 
and ABCDEFGH, is 

Aid 



sinh 



k x V 



a/ 4 At? - Awf ' 

and the hinge ABCD is surrounded by four such boosts, two identical boosts above and 
two below. Thus the final deficit angle measured about ABCD is 



€ 



2 (<>;': ~ <t>?) ■ 



Further details on calculating the deficit angles about a spacelike hinge are 
contained in a recent paper by Brewin [16J. We note that this type of hinge was 
not included in the calculations of Lewis [7], leading to errors in the resulting Regge 
equations. We return to this issue below. 



4. The Regge calculus model 

The vacuum Regge equations take the form p] 
n V" dA * 

where the sum is over all triangles t with area A t that contain the edge Lj, and e t the 

deficit angle about triangle t. Each edge in the lattice yields a single Regge equation, but 

the homogeneity and anisotropy of our model imply that there is one distinct equation 

for each of the four classes of edge in the lattice. 

The Regge equation associated with an individual timelike edge mf involves six 

timelike lattice faces, and has the form 

f)A xt f)A yt r)A zt 

= 2 e xt — + 2 e yt + 2 e zt — 

% dmf 1 dmf % dmf ' 

which is the discrete counterpart of (JIT), the Einstein field equation G tt = 0. Similarly, 
the Regge equation which corresponds to the single spatial edge Ui involves six lattice 
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faces: two spacelike faces with area U{Vi, two with area UiWi, plus the timelike faces 
formed by evolving Ui both forwards and backwards in time. The vacuum Regge 
equation is 

n. _ xy 9A* V dA* z t dAf t dAf_ x 

U " 1 dul + 1 du\ + 1 du\ + du\ 
with similar expressions for the equations arising from the spacelike edges v i and Wi. 

Using the geometric information collected in section [3] we obtain a single Regge 
equation for each class of edge on the ith hypersurface, namely 

Q = (wj + w i+ i) ef _ (vj + Vj+i) ef _ (m + u i+1 ) ef 
4 a/ Aw 2 - 4m 2 ^Avf - 4m 2 4 a/ Aw 2 - 4m 2 



2^/Au^ - 4m?_! 

-2m 2 _ x + ViAv^ 
2 v /A«; i _i-4m?_ 1 

,2 



o = M ^ + + - J- 1 " T 1 ef, + - , ^ ef (9) 
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-2m 2 - 
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- 4mf 


-2m? - 


WjAwj 


2v/A^ 2 


— Amf 



2y/Aw i -.i -4m 2 _ x 

which correspond to the lattice edges m 2 , ttj, Vj and u>j, respectively. 

The structure of the Regge equations (jTjl- tfTOjl is worth considering. Equation (JTJ), 
associated with the timelike edge m 2 , involves edges on, and between, two neighbouring 
hypersurfaces, whereas (18|)-( TT0|) involve information on and between three consecutive 
spatial hypersurfaces. Thus fl7|) is a first-order constraint, while (1H|)-( !T0|) are second-order 
difference equations. 

This contrasts sharply with the equations derived by Lewis [TJ, who neglected the 
curvature associated with spacelike hinges, and was thus unable to derive the spatial 
equations (I8])- (TT0]) . After correctly obtaining (jTJ), Lewis found that he could only make 
sense of the truncated spatial equations by considering a careful average. This averaging 
resulted, once again, in the timelike equation (JTJ. Without the spatial curvature terms 
e( z , ef z and ef y Lewis was unable to build the second-order Regge equations (15])- (fID]) . 

Before examining the continuum limit of the Regge model we consider solutions of 
the discrete equations. Initial data is constructed at to = 1 to match the continuum 
Kasner solution as far as possible. Taking the exact initial data to be x(l) = y(l) = 
z(l) = 1 and x(l) = pi, y(l) = p 2 , and z(l) = p 3 , we mimic the properties of the exact 
solution in the lattice by setting u = v = v = 1 and 

U\ = u + (pi + a) At, vi = v + (p 2 + ol) At, w x = w + (p 3 - 2a) At 

where a is an unknown parameter. This form is chosen to maintain the continuum 
condition x(l) + y(l) + z(l) = 1 to first order, which is physically equivalent to using the 
degrees of freedom in the initial data to generate linear expansion in the volume element. 
With these initial data the Regge constraint (J7J) is solved for the single parameter a. 

A typical solution of the discrete Regge equations is shown in figure [2] for the case 
Pi = 0.75 and At = 0.01. The solution to the initial value problem in this case is 
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(a) Evolution of edge lengths 
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(b) Fractional error in edge lengths 




(c) Time evolution of the constraint Rt (d) Convergence of the mean value of Rt 

Figure 2. Discrete Regge solutions with pi = 0.75 and At = 0.01. 



a = 0.0214379, which represents a roughly 3% change in the initial rate of change of 
Ui compared to the exact estimate. The evolution of the initial data is shown in figure 
while [2b shows the evolution of the fractional error in the Regge solutions compared 
with their exact counterparts. The fractional error in all edges remains in the 5% — 10% 
range, and can be shown to shrink as the time step is reduced. 

The residual error in the first-order Regge constraint ([7]) is defined as 

Rt _ (wj + Wj+^ef | (vi + Vi+i)ef | (uj + u i+1 ) ef 
4 a/ Awf - 4m, 2 4a/A?j 2 - 4m 2 4-y/Aw 2 - 4m 2 ' 

and is a measure of the consistency amongst the Regge equations. Figure [2fc shows 
Rt as a function of time with At = 0.01, and clearly the residual Rt remains small 
throughout the evolution. We repeated this process with different At to estimate the 
rate at which R t reduces as At tends to zero. Figure |2H shows second-order convergence 
in the mean value of Rt over 1 < t < 100 as the timestep is reduced. We explore the 
issue of convergence in more detail in the following sections. 



5. The continuous-time Regge model 



Many of the early applications of Regge calculus to highly symmetric spacetimes 
considered the differential equations that arise in the limit of continuous time and 
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discrete space [TJ El E] • In this section we derive the continuous time Regge equations 
and compare them with the results of Lewis [7] before considering numerical solutions. 

The continuous-time Regge model is developed from the discrete equations in 
section H] in the limit of small At, with the assumption that spatial edges in the lattice 
approach continuous functions of time. For example, the spacelike edge Ui is viewed as 
the value of a continuous function u{t) evaluated at t — t^. A power series expansion 
then relates the edge lengths on neighbouring surfaces, 



du . 1 d 2 u . 9 1 d 3 u . , 
Ui+1 = u{ti) + -At + -—At> + -—At* 



0(At 4 



(12) 



where all derivatives are evaluated at t = U. Similar expressions hold for the edges 
Ui-i,Vi + i etc. The series expansion for the timelike edges m 2 is obtained from ([6]). 
In the continuum limit the deficit angle about the spacelike area A* y is 

^W = J— ^At + 0(At 3 ), 

and the deficit angle about the timelike face formed by the evolution of it, is 



£xt{t) = 27T — 4 cos 



-i 



V w 



+ 0(At), 



^(A-v 2 )(A- w 2 ) 

with similar expressions for the remaining faces. 

Using these expansions, the Regge equations (I7j)- (ITU1) are, to leading order, 
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where e° t , e° t and e° zt denote the zeroth order terms in the deficit angle expansions. The 
leading terms in ffl3|) - ffT6|) are the continuous-time Regge equations, a set of non-linear 
differential equations. Equation f fl3|) is a first-order differential equation, while the 
remainder are second-order. 

The continuous-time Regge differential equations were solved numerically using 
Mathematica. Initial conditions were applied at to = 1> an d the initial data was chosen 
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Figure 3. Solution of the continuous-time Regge differential equations with pi = 0.5. 



to ensure that the expansion of lattice volume elements is initially linear to match the 
exact solution in section EJ Once again introducing a parameter a, we set 

u(l) = i>(l) = w(l) = 1 and 

u(l)=pi + a, v(l)=p 2 + a, w(l)=p 3 — 2a, 

and use the first-order Regge initial value equation (f!3l) to solve for a. This mimics the 
exact initial data, for which a = 0. The second-order, non-linear differential equations 
(|T4l)-(fl6l) are used to evolve the initial data forward in time. 

Figure [3] shows solutions to the continuous-time Regge equations with pi = 0.5. 
The solution to the initial value problem is a = 0.0204161, which represents a small 
deviation (« 4% change in the initial value of u) from the exact initial data. As can 
be seen in the figure, the continuous-time Regge solutions are very similar to the exact 
Einstein solution, with the Regge edges deviating from the exact values by 5% — 7%. In 
the next section we extend the limiting process to the spatial edges, and examine the 
difference between the exact and Regge equations more carefully. 

6. The Regge equations in the limit of continuous space and time 

In this section we explore the discrepancies between the discrete Regge model and the 
Kasner spacetime by examining the truncation error incurred when the Regge equations 
are viewed as approximations of the Kasner-Einstein equations (JX])- (J4j) . We consider the 
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continuum limit of the Regge equations ([7])- (fT0]) as the spatial lattice is refined in both 
space and time. 

The continuous space and time limit of the Regge equations is obtained from the 
temporal series expansions ( TT2|) together with the link between lattice edges and the 
global length scales given by ([5]). Substituting these into the discrete Regge equations 
fl7|)- f|T0|) and simultaneously increasing the number of prisms (n — > oo) while reducing 
the timestep (At — > 0) we obtain a series expansion for the the discrete equations in 
the continuum limit. The refinement parameter n and timestep At are chosen so that 



m 2 < 0. 



The series expansion of the temporal Regge equation ([7]) is 




dy dz dxdz dx dy\ 

x ^i^i + y-n^i + z ^i^i) ( 17 ) 
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^ dz d 2 y ^ dxd 2 z ^ dyd 2 z\ ^ _|_ Q(At 2 ^ ) 
dt dt dt dt dt dt J ' n 2 

which has truncation error that is first-order in the timestep At and second-order in the 
spatial discretization scale 1/n. The spatial Regge equations (1Hi)-(|T0l) are 
At f dy dz d 2 z d 2 y . 9 1 . . 

At (dxdz d 2 z d 2 x 9 1 x 1 

^{*s + ^ + ^ +0(Af V ) } = ° (19) 

At f dx dy d 2 y d 2 x ^, . 9 1.1 „ /n . 
7?UTt +X J + y W +O{At 'n> ) \ =0 (20) 
to leading order in the continuum limit. The truncation error in these equations is 
second order in both the spatial and temporal discretization scales. 

The leading order terms in the continuous time and space Regge equations ( TTTTl - 
(120]) are identical to the Kasner- Einstein equations (EE])- PI, so we expect solutions of 
the discrete lattice equations to approach the continuum solutions as length scales in 
the lattice are reduced. It is clear from the preceding equations that the truncation 
error for the Regge equations, when viewed as approximations to the Kasner-Einstein 
equations, are second order in the spatial discretization scale 1/n. The truncation error 
is also second order in At for the spatial Regge equations (fl8]) - ()2~0]) . 

The truncation error in (IT7j) implies that the Regge initial value equation (J7J) is 
only a first-order approximation to its continuum counterpart. This conflicts with the 
calculations in section H] where the truncation error in the Regge constraint Rt was 
found to converge to zero as the second power of At (see figure [2U). To understand this 
contradiction, we rewrite the coefficients of At in the expansion ( fT7l) as 

3 xyz + y (zx + xz) + z (yx + xy) + x (zy + y'z) 

= x (yz + zy + y'z) + y (xz + z'x + xz) + z (yz + y'x + xy) 

= + O(At 2 ,^), 
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where the final equality follows from substitution of (jl8p - (!20p . Thus the truncation error 
in the Regge constraint f|T7j) is formally of order At, but the coefficient of that term 
in the expansion is a linear combination of the spatial Kasner-Einstein equations. This 
should be zero to leading order for any solution of the spatial Regge equations. 

To clarify this argument, consider again the simulation in section HI Once the 
initial data is set, the numerical evolution is achieved by the repeated solution of the 
spatial Regge equations fl8])- ffT0l) . These are shown above to be second order accurate 
approximations of the Einstein equations in both space and time. We expect from ( jTTI) 
that the leading order error in the Regge constraint is first order in At. However, the 
coefficient of that error term is a linear combination of the Regge equations we are 
solving (to leading order), and thus the coefficient of At is itself zero to second order 
in At. Thus the effective leading order truncation in the Regge constraint equation 
( TT71) is of second order in both space and time. This is consistent with the numerical 
experiments in section HI 

7. Discussion 

In the preceding sections we re-examined one of the most highly symmetric applications 
of Regge calculus to be found in the literature. The primary goal of this study was 
to examine the convergence properties of Regge calculus, and we have shown that for 
the discrete Kasner spacetime the equations of Regge calculus reduce identically to the 
corresponding Einstein equations in the continuum limit. 

The discrete lattice used by Lewis, outlined above, was specifically designed to 
guarantee a one-to-one correspondence between the degrees of freedom in the Regge 
lattice and the metric components in the continuum solution [7]. Despite this, Lewis 
needed to average the Regge evolution equations in order to obtain consistency in the 
continuum limit. The averaging process was chosen to obtain the first order Regge 
equation ([7]), but Lewis was still unable to derive the remaining spatial equations (JED- 
( IT0|) . We showed in section @] that once all lattice curvature elements are included 
in the calculations the Regge equations consist of one constraint and three evolution 
equations. In section [6] we showed that these lattice equations approach the full set of 
Kasner-Einstein equations in the continuum limit. 

It was also shown in section [6] that the discrete Regge equations are second order 
accurate approximations to the Einstein equations for the Kasner cosmology in the limit 
of very fine discretization. This convergence rate is in agreement with many previous 
numerical simulations, in particular the (3+l)-dimensional Regge calculus models of 
the Kasner cosmology that utilized general simplicial lattices [T0| H9] . Unlike the 
current analysis, these simulations did not enforce the homogeneity and anisotropy 
of the Kasner model throughout the evolution, yet displayed second-order convergence 
to the continuum solution. These numerical simulations considered the convergence of 
solutions, rather than equations, and so are complementary to our analysis. 

In general applications of Regge calculus that utilize a simplicial lattice there will 
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be many more Regge equations (one per lattice edge) than Einstein equations (10 per 
spacetime event). The direct comparison between individual Regge and continuum 
equations considered in this paper would not be possible, or even desirable. We expect 
that an appropriate average of the Regge equations would still correspond to the Einstein 
equations in the continuum limit [T7J [TE1 EJ], and several such averaging schemes have 
been suggested. Brewin considered a finite-element integration of the weak-field Einstein 
equations over a simplicial lattice (discussed in [IB]), and suggested that the vertex-based 
equivalent of the vacuum Einstein equations are 

f)A 

j m 

where Ax^ is a vector oriented along edge Lj in a coordinate system based at vertex v. 
The outer summation is over all triangles which meet at v, and the inner sum is over the 
edges on each triangle. These are essentially linear combinations of Regge equations, 
together with boundary terms [T8] . 

Regardless of how one compares the continuum and discrete equations, it is 
ultimately the solutions that are of interest. The application of Regge calculus to the 
Kasner cosmology discussed in this paper demonstrates yet again that Regge calculus is 
a consistent second-order accurate discretization of general relativity, providing further 
support for the use of lattice gravity in numerical relativity and discrete quantum gravity. 
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